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'^ . Abstract 

We present a calculation of the two-loop quark-squark-gluino contributions to pseudoscalar 
Cn ' Higgs boson production via gluon fusion in the MSSM. We regularize the loop integrals using the 

^^ ' Pauli-Villars method, and obtain explicit and compact analytic results based on an expansion in 

the heavy particle masses. Our results - valid when the pseudoscalar Higgs boson is lighter than 
squarks and gluinos - can be easily implemented in computer codes for an efficient and accurate 
determination of the pseudoscalar production cross section. 



e-mail: 

degrassi@f is .uniromaS. it 
divitaOf is .uniromaS. it 
slavichOlpthe .jussieu.fr 



1 Introduction 

With the coming into operation of the Large Hadron Collider (LHC), a new era has begun in the search 
for the Higgs boson(s). At the LHC the main production mechanism for the Standard Model (SM) 
Higgs boson, H^m, is the loop-induced gluon fusion mechanism [1], gg — t- -f^sM; where the coupling 
of the gluons to the Higgs is mediated by loops of colored fermions, primarily the top quark. The 
knowledge of this process in the SM includes the full next-to-leading order (NLO) QCD corrections 
[2, 3, 4], the next-to-next-to- leading order (NNLO) QCD corrections [5] including finite top mass 
effects [6] , soft-gluon resummation effects [7] , an estimate of the next-to-next-to-next-to- leading order 
(NNNLO) QCD effects [8] and also the first-order electroweak corrections [9, 10, 11]. 

The Higgs sector of the Minimal Supersymmetric extension of the Standard Model (MSSM) consists 
of two SU{2) doublets. Hi and H2, whose relative contribution to electroweak symmetry breaking is 
determined by the ratio of vacuum expectation values of their neutral components, tan/3 = vxjv^- The 
spectrum of physical Higgs bosons is richer than in the SM, consisting of two neutral CP-even bosons, 
h and i7, one neutral CP-odd boson, A^ and two charged scalars, H^ . The couplings of the MSSM 
Higgs bosons to matter fermions differ from those of the SM Higgs, and they can be considerably 
enhanced (or suppressed) depending on tan/3. As in the SM case, the gluon-fusion process is one 
of the most important production mechanisms for the neutral Higgs bosons, whose couplings to the 
gluons are mediated by top and bottom quarks and their supersymmetric partners, the stop and 
sbottom squarks. 

In the case of the CP-even bosons h and H the gluon-fusion cross section in the MSSM is known at 
the NLO in QCD.^ The contributions arising from diagrams with quarks and gluons can be obtained 
from the corresponding SM results with an appropriate rescaling of the Higgs-quark couplings. The 
contributions arising from diagrams with squarks and gluons were first computed under the approxi- 
mation of vanishing Higgs mass in ref. [13]. The complete top/stop contributions, including the effects 
of stop mixing and of the two-loop diagrams involving gluinos, were computed under the same approx- 
imation in ref. [14], and the result was cast in a compact analytic form in ref. [15]. Later calculations 
aimed at the inclusion of the full Higgs-mass dependence in the squark-gluon contributions, which are 
now known in a closed analytic form [16, 17, 18, 19]. 

The approximation of vanishing Higgs mass in the contributions of two-loop diagrams allows for 
compact analytic results that can be implemented in computer codes for a fast and efficient evaluation 
of the Higgs production cross section. For what concerns the top-gluon contributions, the effect of 
such approximation on the result for the cross section has been shown [20, 19] to be limited to a few 
percent, as long as the Higgs mass is below the threshold for creation of the massive particles running 
in the diagrams (in this case, the top quarks). While this condition may also apply to the two- loop 
diagrams involving top, stop and gluino, it obviously does not apply to the corresponding diagrams 
involving the bottom quark, whose contribution can be relevant for large values of tan/3. For the 



^First results for the NNLO contributions in the limit of degenerate superparticle masses were presented in ref. [12]. 



latter diagrams the dependence on the Higgs mass should in principle be retained, which has proved 
a rather daunting task. A calculation of the full quark-squark-gluino contributions via a combination 
of analytic and numerical methods was presented in ref. [21] (see also ref. [22]), but neither explicit 
analytic results nor a public computer code have been made available so far. However, ref. [23] 
presented an evaluation of the bottom-sbottom-gluino diagrams based on an asymptotic expansion 
in the large supersymmetric masses that is valid up to and including terms of ©(m^/m^), 0{mb/M) 
and 0(m|/M^), where m^ denotes a Higgs boson mass and M denotes a generic superparticle mass. 
This expansion should provide a good approximation to the full result, at least comparable to the 
one obtained for the top-stop-gluino diagrams, as long as the Higgs boson mass is below all the 
heavy-particle thresholds. An independent calculation of the bottom-sbottom-gluino contributions, 
restricted to the limit of a degenerate superparticle mass spectrum, was also presented in ref. [24], 
confirming the results of ref. [23] . 

In the case of the CP-odd boson A the calculation of the production cross section is somewhat 
less advanced. Due to the structure of the ^-boson coupling to squarks, only loops of top and 
bottom quarks contribute to the cross section at LO, with the bottom loops being dominant for even 
moderately large values of tan/3. In the limit of vanishing yl-boson mass, m^, the contributions from 
diagrams with quarks and gluons were computed at NLO in ref. [25] and at NNLO in ref. [26] (see 
also ref. [27]). For arbitrary values of m^ the NLO contributions arising from two- loop diagrams with 
quarks and gluons, as well as from one-loop diagrams with emission of a real parton, were computed in 
ref. [3]. Supersymmetric particles contribute to the cross section at NLO through two-loop diagrams 
involving quarks, squarks and gluinos. The top-stop-gluino contributions were computed in ref. [28] 
in the limit of vanishing m^. The analytic result for generic values of the stop and gluino masses was 
deemed too voluminous to be explicitly displayed in ref. [28], and was instead made available in the 
fortran code evalcsusy . f [14] . On the other hand, the two-loop bottom-sbottom-gluino contributions, 
which can be relevant for large values of tan/?, have never been directly computed so far. 

In this paper we aim to reduce the gap in accuracy between the available NLO calculations of 
the production cross sections for CP-odd and CP-even Higgs bosons of the MSSM, exploiting the 
techniques we developed for computing the top-stop-gluino [15] and bottom-sbottom-gluino [23] con- 
tributions in the CP-even case. In particular, we present an evaluation of the two-loop top-stop- 
gluino contributions to the pseudoscalar production cross section valid up to and including terms of 
0{m'^/m^) and 0{m'^/M'^). We show how the terms of order zero in m^ can be cast in an extremely 
compact analytic form, fully equivalent to the result of ref. [28], and we investigate the effect of the 
first-order terms. We also evaluate the same contributions via an asymptotic expansion in the large 
superparticle masses, valid up to and including terms of 0{m'j^/M'^) and 0{m^/M'^). While the latter 
result is valid for mt^rriA ^ M but does not assume a hierarchy between mt and m^, the former is 
expected to provide a better approximation in the region with uia < rrit and relatively light superpar- 
ticles, M ~ mt. As a byproduct, we also obtain a result for the bottom-sbottom-gluino contributions 
valid up to and including terms of 0{m'^/m\) and 0{'nn,/M). Finally, we compare our results for the 



bottoni-sbottoni-gluino contributions to both CP-even and CP-odd Higgs production cross sections 
with those obtained in the effective-Lagrangian approximation of refs. [29, 30]. 

A non-trivial technical issue that arises in the calculation of the pseudoscalar production cross 
section is the treatment of the Dirac matrix 75 - an intrinsically four-dimensional object - within 
regularization methods defined in a number of dimensions n^ = 4 — 2e. The original calculation of the 
two-loop quark-gluon contributions of ref. [3] was performed in Dimensional Regularization (DREG), 
employing the 't Hooft-Veltman (HV) prescription [31] for the 75 matrix and introducing a finite 
multiplicative renormalization factor [32] to restore the Ward identities. In ref. [28] the calculation 
of the top-gluon and top-stop-gluino contributions to the Wilson coefficient in the relevant effective 
Lagrangian was performed both in DREG and in Dimensional Reduction (DRED), which, differently 
from DREG, preserves supersymmetry (SUSY). The latter method does not require the introduction 
of finite renormalization factors, but it involves additional subtleties concerning the treatment of the 
Levi-Civita symbol e^^pa- 

In our calculation of the quark-squark-gluino contributions we avoided all problems related to 
the treatment of 75 by employing the Pauli-Villars regularization (PVREG) method. Being defined 
in four dimensions, the PVREG method respects both SUSY and the chiral symmetry, therefore no 
symmetry-restoring renormalization factors need to be introduced. We tested our implementation 
of PVREG by computing the top-gluon contributions via an asymptotic expansion in the top quark 
mass, and recovering the result obtained in DREG in refs. [3, 17]. As a further cross check, we also 
computed the quark-squark-gluino contributions using the DREG procedure outlined in ref. [32], and 
found agreement with the result that we obtained in PVREG. 

The paper is organized as follows: in section 2 we summarize general results on the cross section 
for pseudoscalar Higgs boson production via gluon fusion. In section 3 we outline our implementation 
of the PVREG method. Section 4 contains our explicit results for the NLO contributions arising 
from both top-stop-gluino and bottom-sbottom-gluino diagrams, as well as a discussion of suitable 
renormalization schemes for the bottom contributions and a comparison with the results obtained in 
the effective-Lagrangian approximation. In section 5 we assess the validity of the expansion in powers 
of m^ in the top contributions, and discuss the numerical relevance of the different NLO contributions. 
In the last section we present our conclusions. We also include, for completeness, an Appendix in which 
we present the NLO contributions from one-loop diagrams with emission of a real parton. 



2 Pseudoscalar Higgs boson production via gluon fusion at NLO 

In this section we recall for completeness some general results on pseudoscalar Higgs boson production 
via gluon fusion. The hadronic cross section at center-of-mass energy ^/s can be written as 

o-(/ii + /i2 -^ A + X) = V / dxidx2 faMi^i^f^p) fbMi^'i^f^F) ^ dz5[z —] aab{z) , (1) 

~^Jq Jo \ X1X2J 

where Ta = m\/s, fXp is the factorization scale, fa,hi{x, fXp) the parton density of the colliding hadron 
hi for the parton of type a (for a = g,q,q), and aab the cross section for the partonic subprocess 
ab ^ A + X at the center-of-mass energy s = X1X2S = m^/z. The partonic cross section can be 
written in terms of the LO contribution a^^' and a coefficient function Gab{z) as 

aabiz) = a^''^zGabiz). (2) 

The LO term can be written as 

(0) _ Gf^ag{fiit) 



a(0) = ^^"^^"^ n\' , (3) 

128 ^/2 7r 

where G^ is the muon decay constant and asipn) is the strong gauge coupling expressed in the MS 
renormalization scheme at the scale yU^. Ha is the form factor for the coupling of the pseudoscalar A 
with two gluons, which we decompose in one- and two-loop parts as 

-Ha = n]f + ^nf + 0{al). (4) 

vr 

Due to the structure of the pseudoscalar coupling to squarks (see section 4) , only diagrams involving 
top or bottom quarks contribute to the one- loop form factor Ti]^. The latter can be decomposed into 
top and bottom contributions as 

nl^ = Tf [cot/3/C^^(Tt)+tan/3/C^^(Tb)] , (5) 

where Tp = 1/2 is a color factor, Tq = Arni/m^ , and 



We recall the behavior of K, in the limit in which the pseudoscalar mass is much smaller or much 
larger than twice the mass of the particle running in the loop. In the first case, i.e. r ^ 1, which may 
apply to the top contribution if m^ is relatively small, 

/Ci'(r) -^ -2-^ + 0{t-^), (7) 

while in the opposite case, i.e. r ^ 1, which is relevant for the bottom contribution, 

^"(r) -^ \ ln2(^) + 0{t^) . (8) 
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The analytic continuation of /C(r) corresponds to the replacement mj 



m, 



-ie , thus the imaginary 



part of eq. (8) can be recovered via the replacement ln(— 4/t) — t- ln(4/T) — iir. 

The coefficient function Gab{z) in eq. (2) can be decomposed, up to NLO terms, as 



.(0), 



G'arli^) + ^G^^{^) + Oial), 



(9) 



Gabiz) 
with the LO contribution given only by the gluon-fusion channel: 

Gi?W = S{l-z)6agSbg. (10) 

The NLO terms include, besides the gg channel, also the one-loop induced gq and qq channels: 

G«(z) = 5{l-z) 




+ Pgg{z) In 



/^F 



+ Ca-{1-z + zy Vi{z) + Ca TZgg , 



(11) 



Gqq (z) 



TZr, 



G^'Jiz) 



-qq 1 -- qg 

where the LO Altarelli-Parisi splitting functions are 

1 



Pgqi^) 



ln(l-z) + iln(-| 



+ n, 



qg ' 



PsA^) 



2Ca 



Vo{z) + --2 + z{l-z) 



Pgg{z) = Cf 



1+(1 



(12) 



(13) 



In the equations above Ca = ^c and Cp = {N^ — l)/{2Nc) {N^ being the number of colors), 

/?o = (11 Ca — 2 Nf)/6 (Nj being the number of active flavors) is the one-loop /3-function of the strong 

coupling in the SM, and 

ln*(l-z) 



Vi z) 



1 



(14) 



The two-loop virtual contributions to gg ^ A, regularized by the infrared-singular part of the contri- 
butions from real gluon emission in the one-loop gluon fusion channel, gg — )■ Ag, are displayed in the 
first line of eq. (11). The second line of that equation contains the non-singular contributions from real 
gluon emission. Eq. (12) contains the contributions due to the one-loop quark-antiquark annihilation 
channel, qq — t- Ag, and to the one-loop quark-gluon scattering channel, gq — )■ qA. General expressions 
for the functions TZgg, Tlqq, Tlqg in the case of pseudoscalar production are collected in the Appendix. 
The two-loop form factor ^^^ receives contributions from diagrams involving quarks and gluons, 
as well as from diagrams involving quarks, squarks and gluinos. The contributions from two-loop 
diagrams with quarks and gluons were first computed in ref. [3], and later confirmed in ref. [17]. The 
contribution to li?^ arising from top-stop-gluino diagrams was computed in ref. [28] in the limit of 
vanishing pseudoscalar mass. For what concerns the contribution arising from bottom-sbottom-gluino 
diagrams, no genuine two-loop calculation has been available so far. In the following sections we 
present our calculation of both kinds of quark-squark-gluino contributions. 



3 Technical aspects of the calculation 

In our computation of TiPf we regularized the loop integrals using the PVREG method. For the 
purposes of this computation, the main advantage of PVREG is the fact that all the Lorentz indices 
remain strictly 4-dimensional, thus the 75 matrices anticommute with the other gamma matrices and 
the trace on a string of gamma matrices can be taken using the standard 4-dimensional relations. We 
recall that in PVREG, given an ultraviolet (UV) divergent integral I{q,m?) where q and rn?' denote 
collectively the external momenta and masses, its regularized version is constructed as 

n 

I^iq,m^,Ci,mf) = I{q,m^) + ;^Q/(g,mf). (15) 

In the equation above the original integral I(g, m?) is combined with a number n of replicas, weighted 
by coefficients c^, in which some of the masses of the original integral are replaced by the PV mass 
regulators {rrn), in such a way that the regularized integral is finite if rrii are kept finite, but tends 
to infinity as rrij — t- 00. The number of added terms, as well as the relation that the coefficients Cj 
should satisfy in order to make I convergent, depend on the divergent nature of the original integral. 
If the latter is only logarithmically divergent, a single subtraction is sufficient to construct / , i.e., 
n = 1, ci = —1, mi = Mpy. For what concerns the infrared (IR) divergences associated to massless 
particles, in PVREG they are regularized by giving a fictitious mass A to the massless particle, and 
later considering the limit A — )■ 0. 

All the diagrams contributing to the virtual NLO contributions to pseudoscalar production are at 
most logarithmically UV-divergent, therefore a single subtraction is sufficient to make them convergent. 
In this case, PVREG reduces to subtracting from the original diagrams the same diagrams with 
some of the masses replaced by Mpy, and then taking the limit Mpy — )• 00. In the case of the 
top-gluon contributions also the limit A — )■ must be taken on the fictitious gluon mass. In the 
present calculation, taking the relevant limits for the mass regulators does not introduce additional 
complications with respect to the same calculation performed in DRED or DREG. This is due to the 
fact that we are computing the two-loop diagrams via an asymptotic expansion, so that the final result 
is expressed in terms of two-loop vacuum integrals with different masses and of one-loop integrals. 
Both kinds of terms are fully known analytically, including all the relevant limits when one or more 
masses are sent to infinity or to zero. The asymptotic expansion of the relevant diagrams is generated 
following the procedure described in ref. [23] , which amounts to adding to and subtracting from each 
diagram its IR-divergent part. As discussed in that paper, a diagram minus its IR-divergent part can 
be evaluated via a Taylor expansion in the external momenta (being this combination IR finite by 
construction) while its remaining IR-divergent part, which is expressed as a product of two one-loop 
integrals, must be evaluated exactly. 

In order to test our implementation of PVREG we first considered the two-loop top-gluon contribu- 
tions. These contributions can be split in two parts, one proportional to Cp and the other proportional 
to Ca- The latter, which stems from the non-abelian nature of S'[/(3), is not IR finite but contains a 
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soft and collinear divergence that factorizes on the lowest-order cross-section. In DREG, this IR di- 
vergence appears as a 1/e^ pole multiplying the top contribution to a^^'. We computed the top-gluon 
contributions via an asymptotic expansion in the top mass up to and including terms 0{m^/mf). 
The IR divergences are regularized by giving a mass A to the gluon, while the UV divergences are 
regularized by subtracting to any term a replica in which A is replaced by Mpv- The final result is 
then obtained taking the limits Mpy -^ oo and A — )■ 0. We were able to reproduce in PVREG the 
known result for the top-gluon contributions obtained in DREG [3, 17] once the PVREG IR-divergent 
term 1/2 log (— m^/A^) is identified in DREG with 1/e^. This is quite non-trivial, because it is known 
that, in general, regularizing the IR divergences via a fictitious gluon mass does not respect the non- 
abelian symmetry of 5*^7(3). Thus, one expects to get the correct result only for the part proportional 
to Cp. However, we quantize the Lagrangian employing the Background Field Method (BFM) [33], 
so that the external background gluons satisfy QED-like Ward identities. Then it is not surprising 
that PVREG gives the correct results also for the Ca part. We also remark that within the BFM the 
renormalization of the strong gauge coupling is due only to the wave function renormalization of the 
external background gluons. Thus, the renormalization of ag decouples completely from the rest of 
the calculation, and can be treated separately in the standard way. As a consequence, even if PVREG 
is used to regularize the loop integrals, the LO partonic cross section a^^' can be directly expressed in 
terms of the running coupling as^Hn) as in eq. (3). 

In the evaluation of the top-stop-gluino contributions to "H^ , the two-loop integrals are regularized 
by subtracting from each of them the same expression with m? and m? replaced by Mpy. The 
top-stop-gluino contributions are then computed in two alternative ways: either by means of a Taylor 
expansion in the external momentum, retaining terms of 0{m\/m1) and 0{m\/M'^), or by means of an 
asymptotic expansion in the superparticle masses, retaining terms up to 0{m\/M'^) and 0{m1/M'^). 
The bottom-sbottom-gluino contributions to T-L"^^ can then be recovered from the top-stop-gluino 
contributions computed with the asymptotic expansion, by performing appropriate replacements and 
taking the limit rrn, <C nriA- Considering the hierarchy between mi, and the other masses, we retain 
only terms up to 0{m1/m\) and 0{mij/M). 

We conclude this section with a couple of observations concerning the use of PVREG in the 
computation of the virtual NLO contributions. First, we recall that in PVREG one obtains directly 
the correct result without the need of introducing a finite renormalization factor to restore the Ward 
identities. Second, we note that in PVREG the evaluation of the leading term in the Taylor expansion 
(i.e., the term corresponding to m^ = 0) does not require the computation of counterterm diagrams. 
This seems natural, because the leading term in the one-loop expression, eq. (7), does not depend 
on the top mass. However, the same evaluation in DREG or DRED does require the computation of 
counterterm diagrams. Indeed, in n^ dimensions the one-loop leading term in the Taylor expansion 
contains an 0{e) part that depends on the top mass, so that the counterterm diagrams give rise to a 
non-vanishing contribution. 
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Figure 1: Examples of two- loop quark-gluon diagrams (a), and of two-loop quark-squark-gluino dia- 
grams involving (6) the pseudoscalar-quark coupling or (c) the pseudoscalar-squark coupling. Here, 
q = t,b and i = 1, 2. 



4 Two-loop contributions to the form factor ^^ 

To fix our notation, we write down the Lagrangian for the interactions of the MSSM pseudoscalar A 
with quarks and squarks: ^ 



C D —^htci3 At-f^t -\ ~j^hbSfi Abj^b -\ ^ (^htCi3YtAili2 + hbSf^YbAblb2 
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h.c. 



(16) 



where: ht and /ij, are the top and bottom Yukawa couplings; Yf = A^ — /itan/3 and Y^ = Af, — /xcot/3; 
At and Af, are the soft SUSY-breaking Higgs-squark-squark couplings; fi is the Higgs mass term in 
the MSSM superpotential. Our convention for the sign of ^ is such that, e.g., the stop and sbottom 
left-right mixing angles 9t and Oh obey the relations 
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(17) 



The fact that the pseudoscalar only couples to two different squark mass eigenstates, while gluons 
only couple to two equal eigenstates, implies that the form factor Tij^ receives neither one-loop con- 
tributions from diagrams with squarks nor two-loop contributions from diagrams with squarks and 
gluons. However, contributions to 'H'^ do arise from two-loop diagrams with quarks and gluons, as well 
as from two-loop diagrams with quarks, squarks and gluinos. Examples of such diagrams, involving 
either the pseudoscalar-quark coupling or the pseudoscalar-squark coupling, are given in figure 1. 

The two-loop form factor for pseudoscalar production can be decomposed as 



n 



2£ 



Tf cot p [)Cfg + K 



-ttg 



+ tan/3 {icr, + JClQ 



(18) 



where IC^g denotes the quark-gluon contributions {q = t,b), and /C^^x denotes the quark-squark-gluino 
contributions. In the following we discuss separately the two-loop contributions arising from quark- 
gluon, top-stop-gluino and bottom-sbottom-gluino diagrams. 

^Here and thereafter we use the notation s^ = sin 93, c^ = cos(^ for a generic angle ip. 



4.1 Quark-gluon contributions 



We recall for completeness the results of refs. [3, 17] for the contributions to "H^ arising from diagrams 
with quarks and gluons (see figure la). If the corresponding contribution in the one-loop form factor 



Ti^ is expressed in terms of the physical quark mass, the two-loop contribution for a given quark q 
reads 



/C 



2£ 
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+ CAJ'^{Tq) 



(19) 



If the one-loop form factor is instead expressed in terms of the running quark mass, renormalized in 
the DR scheme at the scale Q, the two-loop contribution becomes 
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(20) 



Expressions for the functions denoted here as J^i(r), J'2{t) and J-3(r), valid for arbitrary values of 



r, can be found in ref. [17]. They correspond to the functions £i ' (4/r) in eq. (4.6), £i ' (4/r 
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in eq. (4.7), and K. 
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(4/r) in eq. (4.12) of that paper, respectively. Their limiting behaviors for 



heavy and light quark are 
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4.2 Top-stop-gluino contributions 



While a fully analytic computation of the top-stop-gluino contributions to 'H.^ valid for arbitrary values 
of all the relevant particle masses is currently beyond our reach, it is possible to derive approximate 
analytic results valid in different phenomenologically relevant limits. 

To start with, we computed the term /C^^- in eq. (18) via a Taylor expansion in the external Higgs 
momentum up to terms of 0{rn\/m^) and ©(m^/M^), where M denotes generically the stop and 



gluino masses. Such expansion should give a reasonable approximation to the full result when tti^ 
is small compared to the other masses, and is anyway restricted to values of m^ below the lowest 
threshold encountered in the diagrams (this usually means ttt-a < 2?7ii). In the limit of vanishing m^ 
we find that our result for /C^^- can be cast in an extremely compact form: 



ttg 



intYt 



rrir 



rrir 



rr 2 2 2 \ St 2 2 2 



(27) 



t2- 



where 



/■/ 2 2 2 



Cf 



+ Ca 



rrit A 



nit 
rrig A 



2/2 2 , 2 \ ^ "H , 2 / 2 , 

m^ [rjig — nit + m^ ) m — ^ + ni^, [rrig + 

9 9/9 9 9 

+ 2 rrig iTit ^{rrig, m^ , m^. 

2/2 2 

nir rrir — m* 

ti ^ ti 1- 

I 2/2, 2 2\^/2 2 2\ 

+ rrig [m^ + m^ — rrig) 'Pyrrig, m^ , m^ ) 



nif 



rrir 



In 



2 

ij 

2 

mi 



2\ 1 nif 2 / 2 

rrix In — tt + rrir rn. 
^ mj *» 



TTlr 

li 



mf 
m?)ln^ 



(28) 



the function (^{m^^ml^rnj ) is given, e.g., in appendix A of ref. [34], and we introduced the shortcut 
A = mf + ni^ + m'j — 2 (m^ m? + m^ m? + m,? m? ) . As appears from eqs. (5) and (7), in the limit of 
vanishing niA the one-loop top contribution to T-La reduces to — cot/3, i.e., it does not actually depend 
on any parameter subject to 0{as) corrections. Therefore, the results in eqs. (27) and (28) do not 
depend on the renormalization scheme in which the calculation is performed. The contributions to 
IC'f^- of the first order in the Taylor expansion in ni^ are too lengthy to be printed here, but in section 
5 we will discuss their relevance in a representative region of the MSSM parameter space. 

The two terms between parentheses in eq. (27) come from the diagrams with pseudoscalar-top and 
pseudoscalar-stop couplings in figures lb and Ic, respectively. Inserting the explicit expressions for 
820^. and Yt we find 



ir2e 

ttg 



mtij. 



rrir 



m 



2- (cot /3 + tan /3) f{m~g,m^,m^J - f{nl~g,m^,m^^) 



(29) 



t2 



i.e., the explicit dependence of /C^|- on At drops out, leaving only a dependence on fi. Ref. [28] 
points out that this happens because the ^ term breaks the axial U(l) Peccei-Quinn symmetry of the 
MSSM potential, thus violating the Adler-Bardeen theorem [35] which would otherwise guarantee the 
cancellation of all contributions from irreducible diagrams beyond one loop. 

We compared our result for /C^/- in the limit of vanishing m^, eqs. (27)-(28), with the result for 

ttg 

the coefficient c{ defined in ref. [28]. That result was deemed too voluminous to be printed explicitly 
in ref. [28], and was made available in the fortran code evalcsusy.f [14]. We find full numerical 



agreement with evalcsusy.f, after taking into account that c| 



;;(!) 



-Tp cot/3/C2| and that ref. [28] 



employs the opposite convention for the sign of ^ with respect to our eq. (17). 
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Even when the superparticles are much heavier than the pseudoscalar, the vahdity of the result 
for /C^^- obtained via a Taylor expansion in tti^ becomes questionable if m-A is close to or even larger 
than rrit. To cover this region of the parameter space we performed an asymptotic expansion of /C^^- 
in the large superparticle masses. More specifically, we consider the case {mA,mt) <C M without 
assuming any hierarchy between ruA and mt, and retain terms up to 0{m\/M'^) and 0{m1/M'^) in 
the expansion. Assuming that the top contribution to T-Ox in eqs. (5) and (6) is expressed in terms of 
the pole top mass, we find 



/C 



2| 
tig 



' 2 ^^'' mt \ 2 m? -m? y \l - xi 



— S2et 1^1 + 



2 m? Yf 



tl t2 



Inxi 



X2 



I- X2 



lnx2 



^t ^ 1 r^U, 



— r^ 2\ '^2 + ^ ^3 - T,^ in) -^ 

rrig (mj — m~ ) ^ ^ ^ 



m 



■^ 7^4 



m| 



7712 — "m' 



2-/-C4 



(30) 



where Xi = mj /ml, the one-loop function K, {t) was defined in eq. (6), and the terms TZi collect 
contributions suppressed by rrit/M or m.^/M'^: 



7^l 



c. 



4(l-xi)3 



+ 



Ca 



2 1 -XI 



:i - X? + 2 xi In xi) ( 2 In ^ - 3 - -/C^^(rt) + 2B 
* mf 2 



>xi Li2(l — xi) — 2xi (3 + xi) Inxi 



:i-xi + xi lnxi)( ln^ + l + i/Ci^(ri)-^ 
w,| 2 



+ 2x1 Li2(l - a:i) + xi (1 + xi) Inxi 



+ 



Cf Yt 



(xi - X2Y m-g 



— l+o^"(Tt 



+ 



Xi 



xf (1-2x2) 
2(1-Xi)(l-X2) ' 2(l-xi)2 



{xi — 2 X2 + xi X2) In xi 



Xl -^ > X2 



(31) 



7^9 



Cf 



4(l-xi) 



2 



2(l-x? + 2xi Inxi) ln3-82:iLi2(l-xi) 



mt 



+ 



Ca 



2(l-xi)2 



+ (1 - x?) ( 1 + llC'Hrt)) - 2x1 ("2 + XI - ^/C1^(tO ) Inxi 



Til 
(1 — Xl + Xl Inxi) In — ^ + 2xi Li2(l — xi) + xi (1 + xi) Inxi 
m| 



Xl i 7- X2 



(32) 



n. 



Cf 



+ 



6(l-xi)4 

Ca 
8(l-xi)3 



;-2-3xi + 6x?-x?-6xi In xi)i 2 + K:^\Tt) -B 
:i-x? + 2xi Inxi)f2 + /C"(ri)-2B) + f xi ^^ X2 



(33) 
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7^4 



Cf Yt 

{xi - X2Y m-g 



+ 



Xl 



xl{l-2x2) 
2(l-xi)(l-a;2) ' 2(1 -xi 



{xi — 2 X2 + xi X2) In xi 



Xl i > X2 



(34) 



In the equations above, B denotes the finite part of the Passarino-Veltnian function BQ^m^ 
computed at the renormahzation scale Q^ = ml 



ml, ml) 
The comparison between the result for /C^^. obtained 



via a Taylor expansion in m^ and the corresponding result obtained via an asymptotic expansion in 
M will be discussed in section 5. 



4.3 Bottom-sbottom-gluino contributions 



A result for the bottom-sbottom-gluino contribution /C^|_ can be obtained by performing the obvious 
replacement t — ;• 6 in the result for /C^^, obtained via the asymptotic expansion in M, eqs. (30)-(34). 
Considering that m^ <^ m,A, and that we are assuming m^ ^ M, we retain only the terms up to 
0{m,b/M) and O^rn^/m.^). In particular, the terms 7^2) ^3 and TZ^ in eq. (30) give contributions of 
higher order in m.b and can be neglected, while in the expression for TZi, eq. (31), we drop the occur- 
rences of K}^{ti,) and use B = 2 — ln(— Tn^/m^). As a result, assuming that the bottom contribution 
to T-L]^ in eqs. (5) and (8) is fully expressed in terms of the pole bottom mass, we again find a rather 
compact expression for the term /C^^_ in eq. (18): 






nibYb 



mf — mi J V 1 — a^i 



Xl . X2 , 

mxi lnx2 



1 - X2 



rrib 



5206^1 • (35) 



Here Xj = m? /mA , and TZi collects the contributions suppressed by rrib/M: 



7^l 



+ 



+ 



Cf 



4(l-xi)3 

Ca 

2(l-xi)2 

Cf Yb 

(xi - X2)^ mg 



(1-Xi -F2xi Inxi) 1- 2 ln( 



-m^ 
ml 



1 — Xl -|- Xl Inxi) ( ln( 



-m^ 
mj 



+ 



Xl 



xf (1-2x2) 

2(1-Xi)(l-X2) ' 2(1-Xi)2 



- 8x1 Li2(l — Xl) — 2xi (3 -|- Xl) Inxi 
1 ) -I- 2xiLi2(l - Xl) -l-xi (1 + xi) Inxi 

Xl — 2 X2 + Xl X2) In Xl 



Xl i > X2 



(36) 



As in the case of the top-stop-gluino contribution, the terms proportional to Yb originate from the 
diagrams that involve the pseudoscalar-sbottom coupling, while the other terms originate from the 
diagrams that involve the pseudoscalar-bottom coupling. Inserting the expressions for 820,^ and Yb in 
the first term in the right-hand side of eq. (35) we obtain 






Cf ^u, ^ m-g^i 

-— -/C (Tb)^ 2 

1 rrir — m,r 

61 b 



(tan /3 + cot /3) 



Xl 



1 — Xl 



Inxi 



X2 



1 - X2 



lnx2 ■S2e, 7^l 



(37) 
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Similarly to what found in ref. [23] for the production of CP-even Higgs bosons, if the one-loop 
contribution to T-L^ is expressed in terms of the pole bottom mass the bottom-sbottom-gluino diagrams 
induce potentially large two-loop contributions. According to whether or not we insert the explicit 
expression for S20[, in our formulae, such contributions manifest themselves either as terms enhanced 
by the ratio nig/mh, as in eq. (35), or as terms enhanced by tan/3, as in eq. (37). However, such 
terms cancel out if the pseudoscalar-bottom coupling entering the one-loop contribution to Tij^ is 
identified with the DR-renormalized mass irif,, while the mass of the bottom quark running in the 
loop is identified with the pole mass Mh (this amounts to rescaling by rhb/Mb the one-loop result fully 
computed in terms of M5). As a result, the two-loop form factor in eq. (18) is shifted as 



n 



2£ 



n 



2i 



\.cm^K.^\n)TpCF 



3 rn| 

4 Q2 



5 ((5mt,) 
4 ra\) 



SUSY 



(38) 



with respect to the result obtained when the one-loop bottom contribution is fully expressed in terms 
of Mfe. Here Q is the scale at which the running mass ifib is renormalized, and ((Jm^)^^^^ denotes the 
SUSY contribution to the bottom self-energy, in units of Cp OLgjix and in the limit of vanishing my, : 



{5mb) 



SUSY 



rub 



Tfl- Tn ~ 



Xl 



1 — Xl 



Inxi 



X2 



X2 



ln2;2 



where 



f{x) 



+ 



x{x — 2) 



In a; . 



(39) 



(40) 



4(l-x) ' 2(1 

While the shift in eq. (38) removes the contributions enhanced by mg/rrib (or tan/3), it does 
introduce potentially large logarithms of the ratio between the renormalization scale Q and the masses 
of the particles running in the loop. Such logarithms cannot be eliminated by a specific scale choice 
for fhii, unless Q is set to a value much smaller than the bottom mass itself. Therefore, as already 
found in ref. [23] for the CP-even Higgs bosons, the bottom contributions to 'H'^ may turn out to 
be sizable even in the "mixed" renormalization scheme in which the tan /3-enhanced contributions are 
absorbed in a redefinition of the pseudoscalar-bottom coupling entering Ti^^. 

Finally, if the bottom contribution to ^^ is fully expressed in terms of the running bottom mass 
fhb the bottom-sbottom-gluino contribution to the form factor in eq. (35) is shifted as 



i-2£ 



3 mi) 



SUSY 



-bbg 



(41) 



In this case 'H'x contains both terms enhanced by rrig/mb and potentially large logarithms, the latter 
arising from (6mb)^^^^ in eq. (41) as well as from the two-loop bottom-gluon contribution in eq. (20). 

4.4 Comparison with the effective-Lagrangian approximation 

It is well known that, in the MSSM, loop diagrams involving superparticles induce interactions be- 
tween the quarks and the "wrong" Higgs doublets, i.e., interactions that are absent from the tree-level 
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Lagrangian due to the requirement that the superpotential be a holomorphic function of the super- 
fields [36]. Such non-holomorphic, loop-induced Higgs-quark interactions result in tan /3- enhanced 
(or tan ;5-suppressed) corrections to the MSSM predictions for various physical observables. If all 
superparticles are considerably heavier than the Higgs bosons they can be integrated out of the La- 
grangian, in which case the loop-induced corrections are resummed in effective Higgs-quark couplings. 
In particular, if (7^ denote the tree-level couplings of a neutral Higgs = {h, H, A) to bottom quarks 
(normalized to the SM value), the corresponding effective couplings gf^ read [29, 30] 



~d 



9i 



l + Afc 



1-A, 



cot a 
tan/3 



~9b 



9^ 



1 + Afe 



1 + A, 



tana 
tan/3 



-at 



9t 



l + Afc 



1 



where a is the mixing angle in the CP-even Higgs sector and, to 0(pi^^ 

as Cp rUg fi tan f3 f xi ^ X2 



A. 



2tt 



m~ 



m~ 



1 — xi 



Inxi 



1 



X2 



\nx2 



Ah cot" /3j , 
(42) 

(43) 



61 62 

In the calculation of processes involving the Higgs-bottom couplings, it is often found that the 
tan /3-enhanced corrections can be included to all orders in an expansion in powers of as tan /3 by 
inserting the effective couplings of eq. (42) in the lowest-order result. A comparison with our explicit 
results for the two-loop form factors allows us to test the validity of that procedure in the case of the 
production of both CP-even [23] and CP-odd Higgs bosons in gluon fusion.^ 



We recall that the bottom-quark contributions Ti 
of the Higgs boson (p = (h, H, A) read 



iifi 



to the one-loop form factors for the production 



n 



U,b 



sma 



cos/3 



7l£ 



^1/2(^6) 



n 



U,b 
H 



Tf 



COS a 



7W 



cos/3 



Si/2in) 



n 



U,b 



U/ 



Tf tannic'' in) , (44) 



,ie,b 



where the function Gi^iT) is given, e.g., in eq. (12) of ref. [23]. Assuming that T-Lj^'" are expressed 
in terms of the pole bottom mass, and that the Higgs-sbottom couplings are renormalized in a way 
that avoids the introduction of additional tan /^-enhanced corrections (see ref. [23] ) , we find that the 
two-loop form factors read 









nj2e 

Ha 



n 



n 



-n 



U,b 
h 

U,b 
H 

U,b 



TT , /, cota\ 


Cf Ab- ficota 2 I X 

+ 4 ^_ S,,^9{X,^X2)^ 

Cf Afo + /i tan a 2 / x 

+ 4 ^_ S,,^9{X„X2)^ 

... , 


+ ... 

+ ... 


, (45) 

, (46) 
(47) 


^b '^ I , n 

as V tan/3/ 

IT , f tana\ 
A , 1 


as \ tan /3 y 

— Afo (1 + cot^ /3) + 
as 



where the ellipses denote contributions suppressed by nih/M or m'^/M'^, as well as all of the contri- 
butions from diagrams involving top and stop, and 



9ixi,X2) 



1 



1 — Xl 



1 + 



Inxi 
1 — Xl 



+ 



1 



1 - X2 



1 + 



Xl 



lnx2 

1 — X2 y Xl — X2 V 1 — 2:1 



Inxi 



X2 



1 - X2 



lnx2 



(48) 



^ A comparison for the light scalar h in the limit of vanishing sbottom mixing was discussed in ref. [24] , and a numerical 
comparison for the heavy scalar H was shown, without a detailed discussion, in ref. [22]. 
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In practice, the effective-Lagrangian approximation consists in rescaling the one-loop bottom con- 
tributions "H , ' by the same factors that rescale the Higgs-bottom coupUngs g'j^ in eq. (42). Expanding 
the rescahng factors to the first order in A5 it is easy to see that the effective-Lagrangian approximation 
does indeed reproduce the two- loop terms proportional to A^ in eqs. (45)-(47). 

It is also interesting to consider the so-called decoupling limit of the MSSM, rriA ^ mz, in which 
cot a —7- — tan /3 and the light scalar h has SM-like couplings to fermions and gauge bosons.^ Eq. (42) 
shows that in this limit the effective coupling of h to bottom quarks is equal to the tree-level coupling, 
therefore in the effective-Lagrangian approximation there are no tan /3-enhanced contributions to Tiff. 
Indeed, for cota — )■ — tan/3 the terms proportional to A^ drop out of the two-loop form factor in 
eq. (45). However, eq. (45) also shows that in the decoupling limit Tiff contains additional tan/3- 
enhanced contributions, controlled by the left-right sbottom mixing X^ = {Ab + fita.n/3), which are 
not reproduced by the effective-Lagrangian approximation. However, when the implicit dependence 
of the sbottom masses and mixing on the bottom mass is taken into account, such contributions turn 
out to be partially suppressed by powers of nib- Indeed, taking for illustrative purposes the limit in 
which the diagonal entries of the sbottom mass matrix as well as the squared gluino mass are all equal 
to M^, and expanding the form factor in powers of nib, we find 

n, D -n^ — -j^ + Tf ^^ + ... , (49) 

where the ellipses denote terms further suppressed by powers of nib or m^, as well as all of the 
contributions from diagrams involving top and stop. The first term in eq. (49) comes from the 
expansion of the terms proportional to s^q in eq. (45), while the second comes from the expansion of 
terms not shown in eq. (45). The contributions neglected by the effective-Lagrangian approximation 
can be relevant for values of Xb large enough to compensate for the suppression due to nib- It should 
however be recalled that in the decoupling limit "H^ ' is not further enhanced by tan/3, therefore - 
differently from what happens in the case of the heavy Higgs bosons - the total form factor for h 
production can still be dominated by the top/stop contributions even for large values of tan/3. 

5 Numerical examples 

We will now illustrate the effect of the two-loop quark-squark-gluino contributions to the form factor 
for pseudoscalar Higgs production in a representative region of the MSSM parameter space. 

The SM parameters entering our calculation include the Z boson mass niz = 91.1876 GeV, the W 
boson mass m^y = 80.399 GeV and the strong coupling constant as{niz) = 0.118 [37]. For the pole 
masses of the top and bottom quarks we take Mt = 173.3 GeV [38] and Mb = 4.49 GeV, the latter 
corresponding to the SM running mass (in the MS scheme) mb{mb) = 4.16 GeV [39]. 



"^The validity of the effective-Lagrangian approximation for the hght scalar h in the decoupling limit was already 
discussed in ref. [30] in the context of Higgs boson decays to bottom quark pairs. 
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Since the squarks do not contribute to the one-loop amphtude for pseudoscalar production, the 
only parameters entering 1-0^ in addition to the quark masses are tan/3 and itia- Neither of those 
parameters is subject to one-loop 0{as) corrections, therefore we need not specify a renormalization 
scheme for them (although it is natural to consider m^ as the pole pseudoscalar mass). The remaining 
input parameters are rrig^ fi, At, Ab and the soft SUSY-breaking mass terms for stop and sbottom 
squarks, niQ, mu and m£). Since these parameters only enter the two-loop part of the form factor we 
need not specify a renormalization scheme for them either. For simplicity, in our numerical examples 
we will set all the SUSY-breaking parameters, as well as the supersymmetric mass parameter /i, to 
a common value M. Note however that the squark mass eigenstates will differ from M, because of 
the supersymmetric (F-term and D-term) contributions to the squark mass matrices as well as of the 
left-right mixing terms. 

In figure 2 we show the top-stop-gluino contribution to the two-loop form factor for pseudoscalar 
production, i.e., the term /C^^- entering eq. (18), as a function of the common SUSY mass M, for 
tua = 150 GeV and tan/3 = 2. Even for the lowest value of M considered in the plot, M = 100 GeV, 
the stop and sbottom masses are above the threshold for real-particle production. The dashed line 
represents the result obtained in the limit of vanishing uia, shown explicitly in eqs. (27) and (28), 
while the solid line represents the result computed at the the first order of the Taylor expansion in 
the pseudoscalar mass, i.e. it includes the effect of terms of 0{m?A/m1) and 0{m\/M'^) which are too 
long to be presented in analytic form. In the computation of these additional terms we assumed that 
the 0{m\/m1) part of the one-loop top contribution, see eq. (7), is expressed in terms of the pole top 
mass. 

It can be seen in figure 2 that the two-loop top-stop-gluino contribution /C^/- is of non-decoupling 

ttg 

nature, i.e., it does not tend to zero when all the superparticle masses become large (note that the 
superpotential parameter /x increases together with the SUSY-breaking parameters). In addition, the 
comparison between the solid and dashed lines shows that when the common SUSY mass M is close 
to TTi-A the combined effect of the terms of 0{m\/m1) and 0{m\/M'^) can be as large as 20%-25% 
with respect to the result obtained for vanishing m^. However, when M increases the effect of the 
terms of 0{m\/M'^) becomes quickly negligible. The remaining discrepancy between the solid and 
dashed lines for moderate to large values of M is due to the terms of ©(m^/m^), and it amounts to 
a modest 6% for the value of m^ considered in this example. 

To assess the importance of the terms of 0{m\/m1) for larger values of m^, we plot in figure 3 
the real part of /C^^- as a function of the pseudoscalar mass, up to a value m^ = 500 GeV well above 
the threshold for real top-quark production. The common SUSY mass is set to the relatively large 
value M = 1 TeV, and tan/3 = 2. As in figure 3, the dashed and solid lines represent the results 
obtained at the zeroth and first order of the Taylor expansion in m\, respectively. The comparison 
between those lines shows that when tha approaches 2 mt the effect of the terms of 0{m\/m^) gets as 
large as 30% with respect to the result obtained for vanishing niA- However, it is natural to wonder 
whether a Taylor expansion in m^ can give an accurate approximation to /C^^- for values of tua close 
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Figure 2: Top-stop-gluino contribution /C^^, as a function of a common SUSY mass Af , for m^ 
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Figure 3: Real part of /Cf^- as a function of m^, for a common SUSY mass Af = 1 TeV and tan/3 = 2. 
The solid and dashed lines are as in figure 2 above, while the dot-dashed line is the result of an 
asymptotic expansion in M which does not assume a specific hierarchy between mt and ttia- 
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to or larger than nif. To address this question, we show in figure 3 as a dot-dashed line the result 
of the asymptotic expansion in M, given explicitly in eqs. (30)-(34). This result was derived under 
the assumption that both m^ and nit are much smaller than M, which is indeed the case for M = 1 
TeV, but it does not require any specific hierarchy between m^ and rrit- The comparison between the 
dot-dashed and solid lines shows that the Taylor expansion at the first order in m^ provides a good 
description of the dependence of /C^^- on the ratio m^/mt up to values of uia of the order of 250 GeV. 
On the other hand, when m^ reaches the threshold for real top production (i.e., at the cusp of the 
dot-dashed line) the result of the asymptotic expansion in M is roughly 80% larger in absolute value 
than the result at the first order of the Taylor expansion in m^, and a full 140% larger than the result 
obtained for vanishing m^. 

In summary, it appears that the compact result for /C^^- given in eqs. (27) and (28), which was 
derived for m^ = 0, can be safely applied only to scenarios in which m^ is smaller than rrit. While 
the inclusion of the terms proportional to m,^ pushes the validity of the Taylor expansion up to larger 
values of niA, the expansion fails when m^ gets close to the threshold for real top production. In 
that case one can use the result of the asymptotic expansion in M, provided that the latter is still 
considerably larger than m^. 

We are now ready to discuss the relative importance of the various two-loop contributions to the 
form factor for pseudoscalar production. We will see that, at least in the region of the parameter space 
that we consider in this example, the results are qualitatively similar to what we found in ref. [23] for 
the case of the heavy scalar H. 

A precise NLO determination of the cross section for pseudoscalar production would require us to 
take into account the contribution of one-loop diagrams with real parton emission, and to perform an 
integration over the phase space (see section 2). However, for the purpose of illustrating the relative 
importance of the various two- loop contributions, we can just define a factor K^ that contains the 
ratio of two-loop to one-loop form factors appearing in eq. (11): 



as „ I'H 



Ka = 1 + 2 ^ Re -f. . (50) 



TT \n 



In the left panel of figure 4 we plot K^ as a function of tan/3, for uia = 150 GeV and all SUSY 
mass parameters equal to M = 500 GeV. The one-loop form factor 1-0^ in eq. (50) contains both the 
top and bottom contributions, computed under the approximations of eqs. (7) and (8), respectively. 
We identify the quark masses in the one-loop form factor with the pole masses, and refer to this 
choice as "on-shell" (OS) scheme. The lines in the plot correspond to different computations of the 
two-loop form factor ^^ : the dotted line includes only the contributions of the top/stop sector (both 
those involving top quarks and gluons and those involving top, stop and gluinos) computed at the 
first order of the Taylor expansion in m\] the dashed line includes also the contributions of two- loop 
diagrams with bottom quarks and gluons; finally, the solid line includes the full contributions of the 
bottom/sbottom sector. 

Comparing the three lines in the left panel of figure 4 it can be seen that the top/stop contributions 
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Figure 4: K factor for the production of a pseudoscalar Higgs yl as a function of tan/3, for mA = 150 
GeV and all SUSY mass parameters equal to M = 500 GeV. The three lines show the effect of the 
different two-loop contributions, in the OS scheme (left panel) and in the "mixed" scheme (right 
panel) . 



dominate the two-loop form factor up to values of tan^ around 5. For larger values of tan/3 the 
contribution of the bottom-sbottom-gluino diagrams (included in the solid line) becomes the dominant 
one, and K^ grows linearly with tan/3. This behavior can be understood by recalling that, as can 
be seen in eq. (16), the Yukawa coupling of the pseudoscalar to bottom quarks is enhanced by tan/3 
with respect to the coupling of the SM Higgs, while the coupling to top quarks is suppressed by tan/3. 
Consequently, for moderate to large values of tan /3 both the one- loop and the two-loop form factors 
in K^ are dominated by the contribution of the diagrams controlled by the pseudoscalar-bottom 
coupling, with the result that the coupling itself cancels out in the ratio. However, the dominant 
contribution from the bottom-sbottom-gluino diagrams in the OS scheme, see eq. (37), contains an 
additional tan/?-enhancement, which explains the linear rise of Kj^. On the other hand, the proximity 
between the dotted and dashed lines shows that, in the OS scheme, the contribution to 'H'j^ of the 
two-loop diagrams with bottom quarks and gluons is very small. This is due to a partial cancellation 
among the three terms entering /C^^ in eq. (19), and to the fact that, in this scheme, the term -F2 (''"&) 
is not enhanced by the potentially large logarithm of the ratio between the bottom mass and the 
renormalization scale, as can be seen by comparing eqs. (19) and (20). 

As discussed in section 4.3, all tan /3-enhanced terms cancel out in a "mixed" renormalization 
scheme in which the pseudoscalar-bottom Yukawa coupling in the one-loop part of the result is identi- 
fied with the DR-renormalized MSSM bottom mass mi,{Q), where Q is a reference scale that we take 
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equal to m^, while the mass of the bottom quark running in the loop is identified with the pole mass 
Ml,. To determine mfe(m^), we first evolve the MS-renormalized SM mass mb{mb) up to the scale 
twa via the NLO-QCD renormalization group equations, then we convert it to the DR-renormalized 
SM mass ml^'^{mjn) via the appropriate shift, and finally we convert it to the MSSM running mass 
according to 

mbiniA) = fhfiniA) , (5f) 

i + Afe 

where A^ is given in eq. (43), and 5b is proportional to the part of {5mbY^^^ in eq. (39) that is not 

enhanced by tan/?: 



4tt 



mj^ m^ — ra^ \1 — Xi l — X2 J 



b-i fe 



(52) 



The "mixed" renormalization prescription is realized by computing the one-loop bottom contribu- 
tion K, {Tb) in eq. (5) in terms of the pole mass Mb, then rescaling it by a factor mb{'mj^)/Mb. The 
two-loop form factor 'H.'^^ must then be shifted as in eq. (38). In the right panel of fig. 4 we present the 
result of this manipulation. The input parameters and the meaning of the different lines are the same 
as for the plot in the left panel. The proximity between the dashed and solid lines, and the flatness 
of the lines for moderate to large values of tan/3, show that the contribution of the two-loop bottom- 
sbottom-gluino diagrams is rather small in this renormalization scheme, and it does not induce an 
additional tan/3-enhancement. However, the comparison between the dotted and dashed lines shows 
that there is a sizable contribution to K^ from the two-loop diagrams involving bottom quarks and 
gluons. This is due to the fact that the shift in eq. (38) brings back a large logarithm, ln(?Ti-^/m^), 
which compensates the scale dependence of the running mass fhb- 

6 Conclusions 

The calculation of the production cross section for the MSSM Higgs bosons is not quite as advanced 
as in the SM. Indeed, despite valiant efforts [21, 22], a full computation of the two-loop quark-squark- 
gluino contributions, valid for arbitrary values of all the relevant particle masses, has not been made 
publicly available so far. Approximate analytic results, however, can be derived if the Higgs bosons 
are somewhat lighter than the squarks and the gluinos. In the MSSM this condition almost certainly 
applies to the lightest scalar h. Moreover, recent results from SUSY searches at the LHC [40] set 
preliminary lower bounds on the squark and gluino masses just below the TeV (albeit for specific 
models of SUSY breaking), suggesting that there might be wide regions of the MSSM parameter space 
in which the condition also applies to the heavy scalar H and to the pseudoscalar A. 

In this paper we presented a calculation of the two-loop quark-squark-gluino contributions to the 
cross section for pseudoscalar production. We exploited techniques developed in our earlier compu- 
tations of the production cross section for the CP-even Higgs bosons of the MSSM [15, 23] to obtain 
explicit and compact analytic results based on expansions in the heavy particle masses. We avoided 
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problems related to the definition of the Dirac matrix 75 in rid ^ 4: dimensions, which are specific to 
the case of pseudoscalar production, by regularizing the loop integrals with the Pauli-Villars method. 
For what concerns the top-stop-gluino contributions, we provided both the result of a Taylor expansion 
in the pseudoscalar mass, up to and including terms of 0{m'j^/mf) and ©(m^/M^), and the result of 
an asymptotic expansion in the superparticle masses, up to and including terms of 0{m\/M'^) and 
0{m'i/M'^). The latter can be easily adapted to the case of the bottom-sbottom-gluino contributions, 
providing a result valid up to and including terms of 0{rn^/m\) and 0{mi,/M). We discussed how the 
tan /^-enhanced terms in the bottom-sbottom-gluino contributions can be eliminated via an appropri- 
ate choice of renormalization scheme for the parameters entering the one-loop part of the calculation, 
and compared our results with those obtained in the effective-Lagrangian approximation. All of our 
results can be easily implemented in computer codes for an efficient and accurate determination of the 
cross section for pseudoscalar production. 

Finally, the results derived in this paper for the production cross section can be straightforwardly 
adapted to the NLO computation of the gluonic and photonic decay widths of the pseudoscalar Higgs 
boson in the MSSM, in analogy to what described in section 5 of ref. [15] for the case of the CP-even 
bosons. 
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Appendix: NLO contributions from real parton emission 

In this appendix we present for completeness our results for the NLO contributions to pseudoscalar 
production from one-loop diagrams with emission of a real parton, i.e., the functions TZgg, TZqq and 
TZqg entering eqs. (11) and (12). Such contributions were first computed in ref. [3] (see also ref. [42]). 
The contribution of the gluon-fusion channel, gg — t- Ag, can be written as 

2 



7^„ 



1 



dv 



"^SSV^' ' ^J 



2\2 



^^^ z{l-z)Jo v{l-v) 
where t = —s{l — z){l — v), u = —s {1 — z) v, and 



m 



W|2 



(l-z + z') 



(Al) 



I'^ggx^j ^' ^/ 



Ti 



cot^ (3 ^** (s, t, u) + tan^ /? AfJs, t, u) 



gg^ 



+ 2 



A%{s,t,u) 



with 



A^qis,t,u) =\A^^is,t,u)\^ + \A'^{u,s,t)\^ + \A^^it,u,s)\^ . 



Defining, for i = t,b , 



m: 



m 



m 



2' 



m: 



2' 



Ui 



2 ' 

mf 



(A2) 
(A3) 

(A4) 



we find: 



\A'^{s,t,u)\' 



ViVj 



6i(s, t, u) H2{si, y.i) Hl{sj,yj) + 62(5, t, u) H2{si,yi) Hl{tj,yj) 

+ h{s, t, u) Hs{si,ti, Ui) Hl{sj,tj,Uj) + 64(5, t, u) H3{si,ti,Ui) hI{uj, Sj,tj) 
+ his, t, u) H2{si,yi) Hl{sj,tj,Uj) + h{s, t, u) H2{si, y,i) Hl{tj,Uj, Sj) 



+ 67(5, t, u) H2{si,yi) Hl{uj,Sj,tj 



+ (i^ J) \ + h.c 



where the function H3{s,t,u) is defined in eq. (2.28) of ref. [19], and 



H2{s,y) = ^ 



log' 



yr^^47^-i 



yi-4/s + iy 

The coefficient functions bi{s,t,u) entering eq. (A5) are 

4*2^,2 



fog' 



^/r^% + i 



(A5) 



(A6) 



bi{s,t,u) 
b2is,t,u) 
63(5, t,u) 
b4{s,t,u) 



(t + u 



2 + s^ - 3tn + s{t + u) + {t + uf 



s^ + t^ + u^ + st + 



2sHu 



2st^u 



{s-t){s + u) {s-t){t + u) 



s^ + t^ + u^ + tu + s{t + u) 
s + t){t + u) , 
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(A7) 

(A8) 

(A9) 

(AlO) 



bQ{s,t,u) 
b7{s,t,u) 



t^ + n^ + s{t + u 
s'^ + {t + u)is + u) + 
s'^ + it + u){s + t) + 



{t — u)ut 



{t + u) 
{u — t)ut 



{t + u) 



The contribution of the quark-antiquark annihilation channel, qq — )■ Ag, can be written as 



n, 



512Z{1-Z) \Aggis 



qq 
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\n 



W|2 
A I 



with 



Aqq{s) = Tf cot 13 ytH2ist,yt) + tan 13 ybH2isb,yb) 



Finally, the contribution of the quark-gluon scattering channel, qg — t- Aq, can be written as 



n, 



qg 



Cf_ ^ /•! dv l + (l-z)V 8^ Aqq{i\ 



i + ji-zY' 

2z 



(All) 
(A12) 
(A13) 

(A14) 
(A15) 

(A16) 



We compared our results for the functions TZgg, TZqq and TZqg with the corresponding results in 
ref. [3], and found full agreement.^ 



^Some misprints in ref. [3] must be taken into account in the comparison. In cq. (C.4) of that paper the term within 
square modulus in the definition of dgq should be divided by 2. Also, the formulae in the Appendices B and C omit all 
occurrences of the MSSM Higgs-quark couplings, denoted in that paper as gg. 
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